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We construct optimal designs for estimating fetal malformation rate, prenatal death rate and 
an overall toxicity index in a toxicology study under a broad range of model assumptions. We 
use WeibuU distributions to model these rates and assume that the number of implants depend 
on the dose level. We study properties of the optimal designs when the intra-litter correlation 
coefficient depends on the dose levels in different ways. Locally optimal designs are found, 
along with robustified versions of the designs that are less sensitive to misspecification in the 
initial values of the model parameters. We also report efficiencies of commonly used designs in 
toxicological experiments and efficiencies of the proposed optimal designs when the true rates 
have non- WeibuU distributions. Optimal design strategies for finding multiple-objective designs 
in toxicology studies are outlined as well. 

Keywords: dose-finding experiment; locally c-optimal design; multiple-objective design; robust 
optimal design; WeibuU model 

1. Introduction 

Developmental toxicity studies play an important role in identifying substances that may 
pose a danger to developing fetuses, including prenatal death and malformation among 
live fetuses. Krewski and Zhu (1995) and Zhu, Krewksi and Ross (1994) demonstrate 
that joint dose-response models for describing prenatal death and fetal malformation in 
developmental toxicity experiments have a good agreement with real data. These models 
can be used to estimate the effective dose corresponding to a gene excess risk for both 
these toxicological end-points, as well as for overall toxicity. It appears that toxicologists 
are generally less receptive to a more rigorous treatment of design issues; see a recent 
foreword/commentary in 2006 in Nature by Giles where the author expounded on the 
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lack of sophistication in current designs for animal experiments. Very recently there 
have been a handful of theoretical articles that utilize statistical principles to design 
toxicology studies. This paper follows this trend and discusses how one may construct 
efficient designs for estimating malformation rate, prenatal death and overall toxicity 
levels under a broad range of model assumptions. A scientifically sound and efficient 
study is crucial because toxicology studies are increasingly more expensive in terms of 
time and labor. An efficient design potentially also means that significantly fewer animals 
will be required in the experiment. In what is to follow, our designs for such experiments 
arc specified in terms of the number of doses to be used, the dose spacing, and the 
proportion of animals to be assigned to each dose. 

Krewski, Smythe and Fung (2002) studied locally optimal designs for the estimating 
the effective dose using joint Weibull dose-response models. The locally optimal designs 
depend on the parameters of the Weibull model and the degree of intra-litter correlation. 
This paper addresses several important design issues in toxicology, such as estimating 
benchmark doses. Estimating benchmark doses has a long history in toxicological studies 
and research continues to this day; some recent papers include (Woutcrscn et al. (2001), 
Moerbeek et al. (2004), Slob et al. (2005)). As in Krewski, Smythe and Fung (2002), 
we seek optimal experimental designs that minimize the variance of the estimated ef- 
fective dose for prenatal death, or malformation or overall toxicity, given the number of 
implants. (Here wc recall an implant is an embryo that has been incorporated into the 
maternal uterus.) These arc two important and common end-points measuring tctrato- 
gcnicity (embryotoxicity) in animal studies (Zhu. Krewski and Ross (1994)). Because 
nonlinear models are employed in such studies, these authors used the simplest design 
strategy and found numerical locally optimal designs for estimating the model models. 
However, it is known that locally optimal designs can depend sensitively on the nominal 
or initial values of the model parameters. In our work, wc present a more sophisticated 
approach to design developmental toxicity experiments so that the implemented design 
is more robust to misspecification in the initial values. Specifically, we first derive analyt- 
ical properties of locally optimal designs for estimating the benchmark dose of prenatal 
death. In particular, we prove several results on the number and levels of doses and in- 
variance properties of the locally optimal designs. We also correct an error in the work 
of Krewski. Smythe and Fung (2002), who used the wrong information matrix for the 
construction of the optimal designs. Second, we study the robustness properties of these 
locally optimal designs with respect to misspecification of the initial parameters. Third, 
we construct locally optimal designs for estimating the effective dose of overall toxicity 
and investigate their sensitivities to misspecification in the initial parameters. Fourth, we 
construct designs that are robust with respect to misspecification of the initial parame- 
ters, and so mitigate a concern raised by some toxicologists. We also investigate relative 
efficiencies of commonly used designs in developmental toxicity experiments. 

Section 2 gives statistical background for our models, which were recently proposed 
in the literature for developmental toxicity studies. In Section 3 we present analytical 
results for locally optimal designs for estimating the effective dose of prenatal death 
and investigate the sensitivity of these designs with respect to misspecification of the 
unknown parameters. Section 4 considers similar problems for estimating the effective 
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dose of overall toxicity, and in Section 5 the methodology is extended to obtain robust 
and efficient designs by a maximin approach. Section 6 evaluates efficiencies of commonly 
used designs in animal studies and briefly discusses efficiencies of optimal designs when 
non-WeibuU probability models are used. All justifications for all our results are deferred 
to the Appendix. 

2. Background for developmental toxicity studies 

In developmental toxicity experiments with laboratory animals such as rats or mice, 
pregnant females are usually exposed to one of several doses of the test agent (including 
a control group at dose zero) during a specified period in gestation. Upon examining 
the uterine contents of each dam, the status of each conceptus is classified and recorded. 
A conceptus may either be dead or alive, and a live fetus may exhibit one or more 
malformations. In industry, the number of litters examined per dose is generally 24 for 
rats, which is one of the two standard species for segment II reproductive studies, at an 
estimated cost of $63,000 for the rat portion of the study. Rabbits are the other standard 
species in a segment II reproductive study and the number of rabbits required per dose 
is usually 12 with the cost of the study of this species being roughly $72,000 (Klaassen 
(2001), page 31-32). 

Let rriij denote the number of implants in the jth litter of the ith dose di, and let 
be the number of prenatal deaths, be the number of live fetuses, and yij be the 
number of fetal malformations. Summary observations from each dam yield a trinomial 
response (r^ , yij , Sy — yij ) conditional on my for which we have 



The fetal malformation rate yij / Sij and the prenatal death rate rij /mij are of particular 
interest. The joint probability of the observed outcome (yij, rij, mij) may be factored as 



where P(my ) is the marginal distribution of the implants number niij. Throughout, we 
let TTi denote the probability of any malformation in a live fetus, 7:2 be the probability 
of the prenatal death and let 0^ be the intra-litter correlation cocfhcicnt within ith dose 
group. Zhu, Krewski and Ross (1994) used generalized estimating equations in conjunc- 
tion with an extended Dirichlet-multinomial covariance function, where the correlation 
coefficient is estimated separately. If = (j/y,ry)"^, the conditional covariance of the 
observation za is 



where fi = 7ri(l — 712), 1/(1 — m^-) < 0^ < 1. 

For simplicity we assume that depends only on the dose level and not on the 
individual litter, that is, rriij =mi= m(di). As pointed out in Krewski, Smythe and Fung 



m^j = Tij + {sij - yij) + yij. 



P{yij,nj,mij) =P(j/y|sy,my)P(ry |my)P(my), 
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(2002) this assumption avoids complicating the model with another level of estimation 
and permits the development of informative designs. Following Zhu, Krewski and Ross 
(1994) we use the Weibull model 

7r,(d) = l-e-"'-'''''^' 

to describe both probabilities tti and 7r2, where d denotes the dose level. Here a^, 6i > 
and 7i > are unknown parameters {i = 1,2). We denote the parameters for the proba- 
bility TTi by 6i. Following Catalano et al. (1993) and Zhu, Krewski and Ross (1994), the 
overall toxicity is defined by 

TT3{d) = l-{l-TT,{d)){l-TT2{d)) (2.1) 

of either a death or malformation occurring. The effective dose EDa for a particular 
probability tt^ is defined as the (unique) solution of the equation 

^(^J^)~^(0) _ 
1-^(0) 

where 7r((i) represents the probability of a response at dose d and a is a given excess risk. 
The excess risk represents additional risk over background among animals that would 
not have responded under control conditions. Typical values for a are 0.01 or 0.05, and 
when a is set at a very low level, say lO"**, the excess risk is also sometimes called the 
benchmark dose or the virtually safe dose (Ryan (1992), Al-Saidy et al. (2003)). Zhu, 
Krewski and Ross (1994) proposed an estimate 9 for estimating 6, the three parameters 
in the Weibull distribution. This estimate is based on quadratic estimating equations 
and has been shown to have reasonable efficiencies for estimating 9. By the (5-method 
(Van der Vaart (1998)) the variance of the estimator EDa for the effective dose can be 
approximated by 

YaiiEDc) « D^Cov{9)D, (2.2) 

and 

D = -^EDa (2.3) 

is the gradient of EDa with respect to 9. We will denote the vector of parameters in tti 
by 9i and its corresponding estimate hy 9i, i = 1,2. 

Throughout, a design is specified by the number of different dose levels, say k, the 
specific dose levels di, . . . ,dk and the proportion of patients, say wi, . . . ,Wk allocated at 
each of these dose levels. In this paper, we consider approximate designs, that is, prob- 
ability measures £, = {rfiiWilfL^ with finite support (Silvey (1980), Pukelsheim (1993)). 
For a given design ^ and total sample size n, the number of observations at each dose 
level Uj is obtained by rounding the quantities nwj to integers, such that X]j=i "-j ~ 
(Pukelsheim and Rieder (1992)). Throughout, we assume for the sake of simplicity that 
the dose range is given by the interval [0, 1] , but the adaption of the methodology to other 
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dose intervals is straightforward. In what is to foUow, we will only present our design 
strategy and results for estimating the effective dose of prenatal death. The strategy for 
estimating the effective dose for a given malformation rate is completely analogous and 
we omit details and corresponding results for space considerations. 



3. Optimal designs for estimating the effective dose of 
prenatal death 



Under the WcibuU model the effective dose for prenatal death equals 

ED„ - 



ln(l - a) ^ ^'"'^ 



b2 

Recalling that 62 = (02:^2:72), the gradient (2.3) in the representation (2.2) is given by 

D-^ED^-~^( l/°2 

"'^^ Vln(-ln(l-a)/fe2)/72, 



{-h-i{l-a)/b2)^/"" 




I/&2 

ln(— ln(l — oi)lh'i)l^2 , 



A = ^7r2(d,) = (l-7r2(d0) I 



72 

li — {di,d2, ■ ■ ■ ,dn]Wi,W2, ■ ■ ■ ,Wn} denotes an approximate design and 

1 

b2dfHd,)^ 

it follows that the covariance matrix of the estimate 62 is approximately 

Cov(§2)-M-\t92), 

where 

" D D'^ 

" §""'"^^(1 + (™. - lmn2{d^)il - 7:2id,)) ^^'^^ 

is the information matrix of the design. Note that the summands in this matrix differ 
by the factors mf from the corresponding terms in the information matrix derived by 
Krewski, Smythe and Fung (2002). There is, in fact, an error in that paper. Consequently 
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we obtain from (2.2) as first-order approximation for the variance of the estimate of the 
effective dose 

Var(fB„) « 02) = D^M-^iC, B^)!) (3.2) 

and a locally optimal design for estimating the effective dose (of prenatal death) mini- 
mizes the function among all designs for which the -E-Dq is estimable. 



It is clear that the information matrix of the optimal design depends on the parameters 
of the model, and, in particular on the quantities nii = in[di) and (pi = (j){di). The simplest 
way to deal with this added complication is to use locally optimal designs proposed 
by Chernoff (1953). This strategy requires that a single prior guess for the unknown 
parameters is available. In developmental toxicity experiments such knowledge is often 
available from preliminary studies. The following results establish properties of locally 
optimal designs for estimating the effective dose of prenatal death. In essence, it says that 
if the excess risk is not too extreme (i.e., near or 1), the locally optimal design requires 
only two doses; otherwise the locally optimal design requires three doses that include the 
extreme levels in the dose interval. The proofs rely on the geometric characterization of 
c-optimal designs of Elfving (1952) and are deferred to the Appendix. 

Theorem 1. Let m and 4> denote the functions defining mi ~ m(di) and (pi = (p{di). If 
the function 

(3 3) 

?7i(d)(l + (m(d)-l)(A(d))7r2(d) ^ ' 

is decreasing, there exist numbers a and a such that the following properties hold: 

(a) If a £ (0,a] U ((5,1), the locally optimal design for estimating the effective dose 
of prenatal death is supported at three points, including the boundary points = and 
d^ = 1 of the design space. 

(b) If a G (a, a), the locally optimal design for estimating the effective dose of prenatal 
death is supported at two points. 

We note that if the function (p is increasing, the assumption of Theorem 1 is satisfied. 
In particular. Bowman, Chen and George (1995) proposed a logistic-type function of the 
form 

Hd) = -, - 1 (3.4) 

for describing the relationship between intra-litter correlation and dose, which is widely 
used in practice. If U2 < this function satisfies the assumptions of Theorem 1. In what 
is to follow, we will use this function to model the intra-class correlation using different 
values of u\ and U2. Of course the case of zero or no correlation is assumed when we set 

U\=U2 = 0. 

The next result tells us that we can limit our search for the optimal design to the 
protocol interval [0, 1] and deduce the corresponding optimal design on the design interval 

[o,r]. 



130 



H. Dette, A. Pepelyshev and W.K. Wong 



Theorem 2. Assume that the quantity rui and the function (f> are constant. The weights 
of the locally optimal design for estimating the effective dose of prenatal death do not 
depend on the parameter 72. Moreover, if d*{a2,b2,"f2,T) are the support points of the 
locally optimal design for estimating the effective dose of prenatal death on the interval 
[0,r], we have 



Our next result shows that if the locaUy optimal design for estimating the prenatal 
death requires three dose levels, these dose levels do not depend on the value of the excess 
risk a. It also provides a complete analytical description of the locally optimal design 
when it is known in advance that the locally optimal design needs only two doses, and 
one of them is the zero dose. 

Theorem 3. Assume that the conditions of Theorem 1 are satisfied. 

(a) The support points of the locally optimal design for estimating the effective dose 
of prenatal death with three support points do not depend on the value of a. 

(b) // the support of a two-point locally optimal design for estimating the effective 
dose of prenatal death contains the point 0, then the second support point is equal to 
EDa and its weight at EDa is equal to W2 ~ g{0)/ {g{0) + g{EDa)) , where 



In Table 1 we display numerical locally optimal designs for estimating the effective dose 
of prenatal death for various combinations of the parameters when the quantity rui and 
the function (f> are assumed to be constants. As stated in Theorem 1 the locally optimal 
designs are either two-point designs or three-point designs because the monotonicity as- 
sumption of the theorem is satisfied. In Table 2 we show some results for a non-constant 
function (j) of the form (3.4), which demonstrate that the assumption of monotonicity on 
the function (3.3) is, in fact, needed. For example, if <j>{d) = 2/(1 + e"^^^) — 1 the corre- 
sponding function in (3.3) is not decreasing. The locally optimal design for estimating 
the effective dose of prenatal death is a 3-point design, but its support dose not contain 
the minimal dose 0. 

In Table 1 and other tables that follow, we also display on the extreme right column 
the efficiency of an equally weighted design supported at five equally spaced points on 
the interval [0, 1]. We denote this design by and note that this is an example of a 
uniform design where doses are equally spread out and an equal number of observations 
are planned at each dose. Uniform designs are intuitively appealing because of their 
simplicity. In practice, approximately uniform designs are used. For example, a recent 
paper by Lee et al. (2006) used 0,0.5,1,2,3,3.5 (mg/kg) dose levels of cadmium, and 
separately, 1,3,5,7.5 dose levels of all-trans-retinoic acid to study their individual and 



d*(a2,fo2,72,l) = (d*(a2,62,l,l))^/^=, 
rft (02 , &2 , 72 , r) = Td* (02 , 62r''" , 72 , 1 ) . 
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combined effects on the induction of forelimb ectrodactyly in C57BL/6 mice. In the cases 
considered in both tables, the results show that this particular uniform design did not 
perform well, averaging about 50%. This means that roughly twice as many rats will be 

Table 1. Locally optimal design for estimating the effective dose of prenatal death conditional 
on the number of implants, assuming the functions and m are constants. The table also 
shows the efficiency of the equidistant design = {0, 1/4, 1/2, 3/4, 1; 1/5, 1/5, 1/5, 1/5, 1/5} (last 
column) 



a 


£12 


b2 


72 


di 


£^2 


£^3 




M2 




ED 




0.05 


0.13 


0.15 


3.33 





0.725 




0.455 


0.545 




0.725 


0.506 


0.05 


0.13 


0.2 


3.33 





0.696 


1 


0.429 


0.545 


0.025 


0.665 


0.539 


0.05 


0.13 


0.25 


3.33 





0.689 


1 


0.404 


0.547 


0.049 


0.621 


0.557 


0.05 


0.13 


0.3 


3.33 





0.682 


1 


0.387 


0.549 


0.064 


0.588 


0.568 


0.05 


0.13 


0.35 


3.33 





0.676 


1 


0.373 


0.551 


0.075 


0.562 


0.575 


0.05 


0.13 


0.4 


3.33 





0.670 


1 


0.363 


0.554 


0.083 


0.540 


0.579 


0.05 


0.01 


0.27 


3.33 





0.607 




0.285 


0.715 




0.607 


0.481 


0.05 


0.05 


0.27 


3.33 





0.653 


1 


0.371 


0.593 


0.036 


0.607 


0.541 


0.05 


0.1 


0.27 


3.33 





0.678 


1 


0.390 


0.558 


0.051 


0.607 


0.558 


0.05 


0.15 


0.27 


3.33 





0.690 


1 


0.399 


0.543 


0.058 


0.607 


0.564 


0.05 


0.2 


0.27 


3.33 





0.698 


1 


0.404 


0.535 


0.061 


0.607 


0.568 


0.05 


0.25 


0.27 


3.33 





0.703 


1 


0.407 


0.529 


0.063 


0.607 


0.570 


0.03 


0.13 


0.27 


3.33 





0.686 


1 


0.367 


0.538 


0.095 


0.519 


0.590 


0.04 


0.13 


0.27 


3.33 





0.686 


1 


0.381 


0.543 


0.076 


0.567 


0.577 


0.05 


0.13 


0.27 


3.33 





0.686 


1 


0.396 


0.548 


0.056 


0.607 


0.562 


0.06 


0.13 


0.27 


3.33 





0.686 


1 


0.412 


0.554 


0.034 


0.642 


0.546 


0.07 


0.13 


0.27 


3.33 





0.686 


1 


0.430 


0.560 


0.010 


0.674 


0.526 


0.08 


0.13 


0.27 


3.33 





0.703 




0.433 


0.567 




0.703 


0.507 


0.1 


0.13 


0.27 


3.33 





0.754 




0.420 


0.580 




0.754 


0.499 



Table 2. Locally optimal design for prenatal death conditional on the number of implants when 
the function m is assumed to be constant and the function <f) modeling the correlation is given 
by (3.4) (a2 — 0.13, i)2 = 0.27, 72 = 3.33, a = 0.05). The table also shows the efficiency of the 
equidistant design ^„ = {0, 1/4, 1/2, 3/4, 1; 1/5, 1/5, 1/5, 1/5, 1/5} (last column) 



Ul 


U2 


di 


£^2 


ds 


Wl 


W2 


W3 


ED 


efr(C") 





-1 





0.636 


1 


0.284 


0.691 


0.025 


0.607 


0.490 





-2 





0.630 


1 


0.242 


0.737 


0.021 


0.607 


0.472 





-3 





0.633 


1 


0.220 


0.757 


0.023 


0.607 


0.464 


-1 


1 


0.082 


0.746 


1 


0.447 


0.482 


0.071 


0.607 


0.588 


-2 


2 


0.071 


0.762 


1 


0.445 


0.487 


0.068 


0.607 


0.569 


-3 


3 


0.052 


0.767 


1 


0.432 


0.503 


0.065 


0.607 


0.551 



132 



H. Dette, A. Pepelyshev and W.K. Wong 



Table 3. Efficiency for estimating tlie effective dose of prenatal deatfi. ^*(^o): locally optimal 
design for 6^ = (a2,''2,72) = (0.13,0.27,3.33) ^„ equidistant design with five different dose lev- 
els (including the largest and smallest dose), and design ^mm = {0, 0.694, 1; 0.349, 0.515, 0.136} 
which is standardized maximin optimal for estimating the effective dose of prenatal death with 
respect to n = [0.05,0.2] x [0.2,0.4] x [2.5,4.5] 



0.2 


0.05 


0.05 


0.05 


0.05 


0.2 


0.2 


0.2 


0.2 


b2 


0.2 


0.2 


0.4 


0.4 


0.2 


0.2 


0.4 


0.4 


72 


2.5 


4.5 


2.5 


4.5 


2.5 


4.5 


2.5 


4.5 




0.802 


0.766 


0.526 


0.967 


0.944 


0.695 


0.749 


0.862 


eff(C.) 


0.522 


0.475 


0.563 


0.521 


0.550 


0.496 


0.588 


0.544 


sff (^mm ) 


0.808 


0.740 


0.663 


0.923 


0.920 


0.665 


0.872 


0.822 



needed in the uniform design to obtain estimates for the parameters as accurate as those 
provided by the locally optimal design. 

In general, our numerical results show that there are four types of locally optimal 
designs for estimating the effective dose of prenatal death, namely: 

{0,d2,l;wi,W2,W3}, {0,(^2; wi,W2}, {di,d2, 1; wi,u'2, W3}, {di,d2;wi,W2}. 

Moreover, if the assumptions of Theorem 1 are satisfied, there exist only two types, that 
is, 

{0,d2,l;wi,W2,W3}, {di,d2;wi,W2}. 

Before any design is implemented, it is useful to investigate the robustness of the 
locally optimal designs for estimating the effective dose with respect to misspecifica- 
tion in the initial parameters. For this purpose we consider the locally optimal $,*{Oo) = 
{0,0.686,1; 0.396, 0.548, 0.056} for the parameter 9'^ = (02,62,72) = (0.13,0.27,3.33) and 
calculate the efficiency 

em) = (3.5) 

for various values of 6. These results arc listed in Table 3. We observe that lo- 
cally optimal designs are not too sensitive with respect to changes of the param- 
eter a2, but a misspecification of the parameters 62 and 72 has more serious ef- 
fects. The table also shows the corresponding efficiencies of the equidistant design 
= {0, 1/4, 1/2, 3/4, 1; 1/5, 1/5, 1/5, 1/5, 1/5}. In most cases these are smaller than the 
efficiencies of the locally optimal design for estimating the effective dose. In addition, the 
table contains efficiencies of a maximin design Cmm, whose construction will be motivated 
in Section 5. This design performs substantially better than the uniform design ^„ and 
achieves nearly the same efficiencies as the locally optimal design £,*{0o) in those cases 
where C*(6'o) is very efficient. 
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4. Dose-finding for overall toxicity conditional 
number of implants 

If two WcibuU models with parameters Of = (ai, 61, 71) and 62 = {02, 62, 72) arc used for 
modeling the overall toxicity in (2.1), the effective dose based on ^^{d) is defined as a 
solution of the equation 



or, equivalently. 



a = 1 - exp{biEDl' + b2EDZ^), 



\n{l-a) = biEDl' +b2EDl\ 



The approximation for the variance of the estimator based on generalized estimating 
equations is given by (2.2), where 9'^ ~ (^1,^2) and 



- d -1 

de bi-fiEDZ'-' + b2j2EDZ'-^ 



( \ 

EDl^ 

biEDl^ \n{ED^) 


EDl^ 
\b2EDlHn{EDo,)J 



If ^ = {di,d2, . . . ,dn','Wi,W2, ■ ■ ■ jWn} denotcs an approximate design we have 

Cov(^)«A/-i(e,0), 
where the information matrix is given by 

and the two non-vanishing blocks are defined by 



Mi(e,^) = E 



^ Ya.T{y^\m^) 



E 



M2itO)=J2 



^ 'm^{l + (m, - l)0O7ri(d,)(l -^2(dO)(l - Mdi){l - Md^)))' 

^(2)j^f2)» 



^ Var(ri|TOi) 



E 



^ TOi(l + (mj - l)4>i)n2idi){l - n2{d^)) ' 
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with 

DU). = 4^^3{d^) - (1 - ^^Ad^)) I dl^ I , J 1,2. 

Note that -M(^, 6') is a block-diagonal matrix and, as a consequence, the optimality crite- 
rion minimizing the variance of the estimate for EDa can be interpreted as a composite 
optimality criterion in the sense of Lauter (1974), that is, 

= bl^M^\i,6)b(,)+bl^M^\i,6)b(^y (4.1) 

It is intuitively clear that locally optimal designs for estimating the EDa for overall 
toxicity are three-point designs if the parameters in 'ni{d) and ■n2{d) are similar. In all 
cases of practical interest these designs have to be calculated numerically. Some examples 
of optimal designs are presented in Table 4 for constant functions rrii and (pi. Table 5 
presents optimal designs for the case where the correlation is of the form (3.4) and 
the mi's are constants. We observe that in most cases, the locally optimal designs are 
supported at three points, but there are also situations (in particular for large differences 
between the parameters 71 and 72), where four different dose levels arc required for 
the optimal estimation of the effective dose of the overall toxicity. The results of our 
investigation of the robustness properties of the locally optimal designs for estimating the 
effective dose with respect to misspecification of the initial parameters are summarized 
in Table 6. 

We next investigate whether the locally optimal design for estimating the effective dose 
of prenatal death is efficient for estimating the effective dose of overall toxicity. We also 
compare the optimal design with the equidistant design with five different dose levels. 
In Tables 7 and 8, we display efficiencies of the two designs for various combinations 
of the parameter 9 to study their robustness for estimating the effective dose of overall 
toxicity when the initial parameters have been misspecified and the design is optimal for 
estimating the effective dose of prenatal death. 

We observe that the performance of the locally optimal design for estimating the 
effective dose of overall toxicity depends sensitively on changes of the parameters 61 
and 62. If bi is very different from the parameter 62 used in the construction of the 
locally optimal design for estimating the effective dose of prenatal death, this design 
becomes inefficient for estimating overall toxicity. In such cases even the uniform design 
performs better. Otherwise the locally optimal design for estimating the effective dose of 
prenatal death is at least as good as the uniform design (and in many cases substantially 
better). The table also shows efffciencies of the design ^mm, which will be constructed 
in the following section as a robust and efficient alternative to locally optimal designs. 
The design S^mm performs uniformly better than the locally optimal design ^*(0o) for 
estimating the effective dose of prenatal death. In many cases, it is substantially more 
efficient than the uniform design and in the cases where the equal allocation rule yields 
the best efficiencies, the loss of efficiency obtained from S^mm is rather small. 
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Table 4. Locally optimal designs for estimating the effective dose of overall toxicity conditional 
on the number of implants. The functions m and (f) are constant, a — 0.05 and denotes the 
equidistant design with five different dose levels 0, 1/4, 1/2,3/4, 1 



ai 


bi 


71 


a2 


62 


72 


di 


rf2 


d3 


d4 


Wl 


W2 


W3 


U14 


eff(C«) 


0.06 


0.7 


2 


0.13 


0.15 


2 





0.495 


1 




0.330 


0.546 


0.124 




0.653 


0.06 


0.7 


2 


0.13 


0.15 


3.33 





0.493 


1 




0.291 


0.574 


0.134 




0.699 


0.06 


0.7 


3.37 


0.13 


0.15 


2 





0.573 


1 




0.372 


0.536 


0.092 




0.593 


0.06 


0.7 


3.37 


0.13 


0.15 


3.33 





0.658 


1 




0.331 


0.546 


0.123 




0.634 


0.06 


0.5 


3.37 


0.13 


0.3 


3.33 





0.665 


1 




0.333 


0.549 


0.118 




0.607 


0.06 


0.7 


3.37 


0.13 


0.3 


3.33 





0.653 


1 




0.321 


0.551 


0.128 




0.619 


0.06 


0.9 


3.37 


0.13 


0.3 


3.33 





0.640 


1 




0.311 


0.551 


0.138 




0.630 


0.06 


0.7 


3.37 


0.05 


0.3 


3.33 





0.630 


1 




0.299 


0.577 


0.124 




0.593 


0.06 


0.7 


3.37 


0.25 


0.3 


3.33 





0.673 


1 




0.338 


0.532 


0.130 




0.635 


0.02 


0.7 


3.37 


0.13 


0.3 


3.33 





0.646 


1 




0.315 


0.559 


0.126 




0.641 


0.09 


0.7 


3.37 


0.13 


0.3 


3.33 





0.657 


1 




0.325 


0.546 


0.129 




0.613 


0.02 


1.2 


2.2 


0.05 


0.2 


3.7 





0.402 


0.636 


1 


0.212 


0.620 


0.040 


0.129 


0.655 


0.02 


1.2 


2.2 


0.05 


0.2 


3.3 





0.421 


0.541 


1 


0.221 


0.596 


0.041 


0.141 


0.680 


0.02 


0.9 


2.2 


0.05 


0.2 


3.7 





0.434 


0.590 


1 


0.228 


0.585 


0.065 


0.121 


0.674 


0.02 


1.6 


2.2 


0.05 


0.2 


3.7 





0.368 


0.713 


1 


0.198 


0.636 


0.029 


0.136 


0.632 



Table 5. Locally optimal designs for estimating the effective dose of overall toxicity conditional 
on the number of implants. The functions m is constant, while the correlation function (j) is given 
by (3.4), ai = 0.06, 61 = 0.7, 71 = 3.37, a2 = 0.13, 62 = 0.3, 72 = 3.33, a = 0.05 and C« denotes 
the equidistant design with five different dose levels 0, 1/4, 1/2,3/4, 1 



Ul 


U2 


di 


d2 


d3 


Wl 


1^2 


W3 







-1 





0.600 


1 


0.227 


0.652 


0.121 


0.554 





-2 





0.594 


1 


0.192 


0.690 


0.118 


0.538 





-3 





0.596 


1 


0.174 


0.708 


0.118 


0.530 



5. Robust and efficient designs for prenatal death 

As pointed out in the previous sections, locally optimal designs are not necessarily robust 
with respect to a misspecification of the unknown parameters. To obtain designs that 
are efficient and robust over a certain range of the parameters for the WeibuU model, we 
study a maximin approach proposed by Miiller (1995) and Dette (1997), which assumes 
that there is prior information on the range of plausible values of unknown parameters. 
To be precise, we concentrate on optimal designs for estimating the effective dose of 
prenatal death, where the correlation function is given by the one parametric logistic 
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family 

^W = T^^-i- (5-1) 

We assume that the experimenter has some knowledge about the location of the param- 
eters, that is, 

02 e [a,a], 62efe6], 72 e [TiT]; 

For given 6^ = (02,^2,72) and u, define £,*{9,u) as the locally optimal designs for esti- 
mating the effective dose and, for a given design, define 

""^^^ ' DTM-^{^,e,u)D ^ ' 

A design is called standardized maximin optimal for estimating the effective dose if 
it maximizes the worst efficiency over some set of the parameters, that is, 

^mm = argmax min efl:g^ii(^, 6', u). (5.3) 

J {e,M)£0 

Here the set f2 is defined by $7 = [a, a] x [6, h] x [7, 7] x [u, m] and is user-selected. Optimal 
designs with respect to this robust criterion have to be calculated numerically in all cases 
of practical interest. In Table 9 we display standardized maximin optimal designs with 
respect to various sets VL assuming the quantities and the correlation function (5.1) 
are constants, that is, u = u = 0. We observe that in all situations the standardized 
maximin optimal designs are supported at three points and they include the largest and 
smallest doses. However, the results of Braess and Dette (2007) indicate that there will 
also exist standardized maximin optimal designs for estimating the effective dose with 



Table 6. Efficiency for estimating the effective dose of overall toxicity using 
three designs: ^*(So), the locally optimal design for 9q = (ai, fei, 71, 02, &2, 72)"^ = 
(0.06, 0.7, 3.37, 0.13, 0.27, 3.33), the equidistant design with five different dose levels (including 
the largest and smallest dose) and the design S^mm ~ {0, 0.694, 1; 0.349, 0.515, 0.136}, which is 
standardized maximin optimal for estimating prenatal death with [0.05, 0.2] x [0.2, 0.4] x [2.5, 4.5] 



ai 


0.03 


0.03 


0.03 


0.03 


0.09 


0.09 


0.09 


0.09 


bi 


0.4 


0.4 


0.9 


0.9 


0.4 


0.4 


0.9 


0.9 


71 


2.7 


4.2 


2.7 


4.2 


2.7 


4.2 


2.7 


4.2 


a2 


0.05 


0.05 


0.05 


0.05 


0.2 


0.2 


0.2 


0.2 


b2 


0.2 


0.2 


0.4 


0.4 


0.2 


0.2 


0.4 


0.4 


72 


2.5 


4.5 


2.5 


4.5 


2.5 


4.5 


2.5 


4.5 


eff(r(eo)) 


0.875 


0.873 


0.681 


0.982 


0.981 


0.705 


0.880 


0.882 


eff(Cn) 


0.588 


0.566 


0.594 


0.584 


0.606 


0.574 


0.619 


0.614 


sff (^mm ) 


0.730 


0.965 


0.531 


0.942 


0.906 


0.871 


0.743 


0.984 
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Table 7. Efficiency for estimating the effective dose of overall toxicity. ^*{0o) = 
{0, 0.686, 1; 0.396, 0.548, 0.056}: locally optimal design for estimating the effective dose of pre- 
natal death {Gq = (02,62,72)"^ = (0.13,0.27,3.33), constant correlation), equidistant design 
with five different dose levels 0, 1/4, 1/2, 3/4, 1 and design ^mm = {0, 0.694, 1; 0.349, 0.515, 0.136}, 
which is standardized maximin optimal for estimating the effective dose of prenatal death with 
respect to « = [0.05,0.2] x [0.2,0.4] x [2.5,4.5] 



ai 


0.02 


0.02 


0.02 


0.02 


0.1 


0.1 


0.1 


0.1 


bi 


0.3 


0.3 


1.1 


1.1 


0.3 


0.3 


1.1 


1.1 


71 


2.5 


4.5 


2.5 


4.5 


2.5 


4.5 


2.5 


4.5 


eff(r(eo)) 


0.762 


0.977 


0.269 


0.883 


0.815 


0.977 


0.328 


0.899 


eff(C.) 


0.669 


0.583 


0.770 


0.611 


0.618 


0.594 


0.644 


0.598 


sff (^7nm ) 


0.901 


0.986 


0.437 


0.984 


0.935 


0.980 


0.520 


0.995 


ai 


0.03 


0.03 


0.03 


0.03 


0.09 


0.09 


0.09 


0.09 


bi 


0.4 


0.4 


0.9 


0.9 


0.4 


0.4 


0.9 


0.9 


71 


2.7 


4.2 


2.7 


4.2 


2.7 


4.2 


2.7 


4.2 


eff(r(eo)) 


0.754 


0.959 


0.443 


0.885 


0.796 


0.965 


0.491 


0.899 


eff(Cn) 


0.652 


0.589 


0.708 


0.609 


0.622 


0.592 


0.646 


0.599 


sff (^mm ) 


0.902 


0.991 


0.646 


0.984 


0.931 


0.991 


0.703 


0.993 



a larger number of support points. The table also contains the minimal efficiency of the 
standardized maximin optimal design for estimating the effective dose, that is, 

min effau(e,6',u). 



Table 8. Efficiency for estimating the effective dose of overall toxicity. 5* (60) ~ 
{0, 0.686, 1; 0.396, 0.548, 0.056}: locally optimal design for estimating the effective dose of pre- 
natal death {6q = (02,62,72)"^ = (0.13,0.27,3.33), constant correlation), equidistant design 
with five diflferent dose levels 0, 1/4, 1/2, 3/4, 1 and design = {0, 0.694, 1; 0.349, 0.515, 0.136} 
which is standardized maximin optimal for estimating the effective dose of prenatal death with 
respect to n= [0.05,0.2] x [0.2,0.4] x [2.5,4.5] 



ai 


0.03 


0.03 


0.03 


0.03 


0.09 


0.09 


0.09 


0.09 


bi 


0.4 


0.4 


0.9 


0.9 


0.4 


0.4 


0.9 


0.9 


71 


2.7 


4.2 


2.7 


4.2 


2.7 


4.2 


2.7 


4.2 


a2 


0.05 


0.05 


0.05 


0.05 


0.2 


0.2 


0.2 


0.2 


b2 


0.2 


0.2 


0.4 


0.4 


0.2 


0.2 


0.4 


0.4 


72 


2.5 


4.5 


2.5 


4.5 


2.5 


4.5 


2.5 


4.5 


eff(r(eo)) 


0.570 


0.968 


0.355 


0.821 


0.767 


0.892 


0.532 


0.914 


eff(Cn) 


0.588 


0.566 


0.594 


0.584 


0.606 


0.574 


0.619 


0.614 


sff (^mm ) 


0.730 


0.965 


0.531 


0.942 


0.906 


0.871 


0.743 


0.984 



138 



H. Dette, A. Pepelyshev and W.K. Wong 



Table 9. Standardized maximin optimal designs for estimating the effective dose of prenatal 
death conditional on the number of implants. The functions m and 4> are constant, a = 0.05 and 
^„ denotes the equidistant design with five different dose levels 0, 1/4, 1/2, 3/4, 1 



9l2 


02 




62 


I2 


72 


di 


d2 


d3 




W2 


W3 


min off 


min off 


0.1 


0.12 


0.25 


0.3 


3.1 


3.5 





0.681 


1 


0.387 


0.551 


0.063 


0.976 


0.549 


0.1 


0.15 


0.25 


0.3 


3.1 


3.5 





0.684 


1 


0.389 


0.546 


0.065 


0.972 


0.549 


0.1 


0.17 


0.22 


0.3 


3 


3.7 





0.696 


1 


0.390 


0.536 


0.073 


0.930 


0.533 


0.08 


0.18 


0.21 


0.33 


2.6 


4 





0.691 


1 


0.371 


0.524 


0.105 


0.809 


0.514 


0.07 


0.19 


0.2 


0.34 


2.5 


4.1 





0.690 


1 


0.366 


0.519 


0.115 


0.757 


0.502 



Table 10. Standardized maximin efficient optimal design for prenatal death conditional on the 
number of implants. The function m is constant, while the correlation function is given by (5.1) 
and a = 0.05. denotes the equidistant design with five different dose levels 0, 1/4, 1/2, 3/4, 1 
and various sets are considered in the standardized maximin optimality criterion (5.3), Q.\{u,u) = 



[0.07,0.19] 


X 


[0.19,0.34] X [2.5,4.1] x 


[u. It 


], n2{u, 


«) = [( 


0.1,0.12] 


X [0.25,0.3] x [3.1,3.5] 


X [u,u] 




u 


u 


di 


d2 


da 


di 


Wl 


'W2 


W3 


W4, 


mineff(^*) mineff(^„) 







1 





0.469 


0.721 


1 


0.269 


0.258 


0.400 


0.073 


0.654 


0.412 







2 





0.460 


0.722 


1 


0.232 


0.282 


0.417 


0.069 


0.640 


0.392 




1 


2 





0.545 


0.665 


1 


0.239 


0.110 


0.535 


0.117 


0.655 


0.392 







1 





0.653 


1 




0.343 


0.599 


0.058 




0.896 


0.469 


0,2 {u, u) 





2 





0.649 


1 




0.327 


0.615 


0.057 




0.873 


0.449 




1 


2 





0.637 


1 




0.254 


0.700 


0.046 




0.946 


0.449 



and the minimal efficiencies of an equidistant design with five different dose levels. Note 
that the standardized maximin optimal designs yield reasonable efficiencies over the full 
set n and that the minimal efficiency of the uniform design over this set is substan- 
tially smaller. In Table 10 we consider the case where the correlation can be modeled by 
the function (5.1). Wc observe that the standardized maximin optimal designs are sup- 
ported at three or four points and, compared to Table 9, the efficiencies are smaller. This 
is intuitively clear because wc have incorporated more robustness with respect to the 
assumption of a constant correlation in the construction of efficient designs for estimat- 
ing the effective dose. Again the equidistant design yields substantially smaller minimal 
efficiencies compared to the standardized maximin optimal design. 

6. Efficiency of standard designs and concluding 
remarks 

It is interesting to evaluate the efficiencies of commonly used designs in developmental 
toxicity studies. One such class is the set of uniform designs. These designs are equally 
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spread out in the dose interval of interest and an equal number of animals is allocated 
to each dose. As such, they are intuitive and easy to implement. Krewski, Smythc and 
Fung (2002) provided an overview of experimental designs for 11 developmental toxicity 
studies conducted under the U.S. National Toxicology Program. In their Tabic 1, they 
listed the doses employed in these studies that involved either rabbits, rats or mice. The 
designs usually have a roughly equal number of animals at each dose and some of their 
dose levels, after scaling to our protocol interval [0,1], are listed in our Tables 11 and 12. 
Following Krewski. Smythc and Fung, we call these "standard" designs. 

Tables 11 and 12 display the efRciencies of "standard" designs for estimating the pre- 
natal death rates and the overall toxicity rate. We observe that the standard design can 
perform poorly when model parameters are misspecified. For instance, the efficiencies 
of the standard design listed in the first row can be less than 30% for estimating the 
prenatal death rate and the overall toxicity rate. Some standard designs have efRciencies 
as low as 0.22 for estimating the prenatal death rate. Interestingly, the uniform design 
with five doses has at least 50% for all cases shown in the tables. The second parts of the 
two tables show the efficiencies of a uniform design with three support points (note that 
the local optimal design usually has two or three support points). One observes that this 
design yields very low efRciencies and cannot be recommended. In general, it is advisable 
that the researcher assess the efRciencies of a design under different optimality criteria 
before its implementation. 

In practice, there are usually several objectives in the study and these objectives may 
not be of equal interest to the researcher. For instance, the researcher may be interested 
in designing a study whose primary aim is to estimate the prenatal death rate, secondary 



Table 11. Efficiency of "standard" designs for estimating EDa for prenatal death with different 
values of parameters, <j)(d) = constant 













a2 


0.13 


0.05 


0.13 


0.13 


0.05 


0.13 














0.27 


0.27 


0.15 


0.27 


0.15 


0.27 




d2 


ds 


d4 




72 


3.3 


3.3 


3.3 


2 




2 


1 





0.25 


0.5 


1 






0.28 


0.31 


0.23 


0, 


.53 


0.46 


0.75 





0.33 


0.67 


0.83 


1 




0.61 


0.55 


0.57 


0, 


.53 


0.48 


0.52 





0.25 


0.5 


0.75 


1 




0.56 


0.54 


0.51 


0, 


.57 


0.51 


0.62 





0.3 


0.5 


0.7 


1 




0.54 


0.54 


0.46 


0, 


.62 


0.54 


0.65 





0.17 


0.33 


0.67 


1 




0.49 


0.47 


0.44 


0, 


.50 


0.45 


0.63 





0.05 


0.15 


0.5 


1 




0.28 


0.31 


0.24 


0, 


.49 


0.43 


0.51 





0.125 


0.25 


0.5 


1 




0.26 


0.29 


0.22 


0, 


.46 


0.40 


0.64 





0.1 


0.2 


0.5 


1 




0.27 


0.30 


0.23 


0, 


.46 


0.40 


0.58 





0.3 


1 








0.04 


0.05 


0.03 


0, 


.29 


0.27 


0.73 





0.4 


1 








0.14 


0.18 


0.11 


0, 


.52 


0.45 


0.71 





0.5 


1 








0.34 


0.39 


0.27 


0, 


.69 


0.60 


0.58 





0.6 


1 








0.58 


0.61 


0.47 


0, 


.73 


0.64 


0.40 





0.7 


1 








0.74 


0.69 


0.64 


0, 


.59 


0.53 


0.23 
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Table 12. Efficiency of "standard" designs for estimating EDa for overall toxicity with different 
values of parameters in the WeibuU model with (j){d) = constant 













fli 


0.06 


0.06 


0.06 


0.06 


0.06 


0.06 


0.02 












bi 


0.7 


0.7 


0.7 


0.7 


0.7 


0.2 


0.7 












71 


3.37 


3.37 


3.37 


3.37 


1 


3.37 


3.37 












£12 


0.13 


0.05 


0.13 


0.13 


0.13 


0.13 


0.13 












b2 


0.3 


0.3 


0.1 


0.3 


0.3 


0.3 


0.3 


di 


d2 


d-i 


d4 


d5 


72 


3.33 


3.33 


3.33 


1 




3.33 


3.33 


3.33 





0.25 


0.5 


1 






0.36 


0.39 


0.34 


0, 


,78 


0.69 


0.29 


0.37 





0.33 


0.67 


0.83 


1 




0.61 


0.56 


0.64 


0, 


.53 


0.40 


0.63 


0.63 





0.25 


0.5 


0.75 


1 




0.62 


0.59 


0.64 


0, 


.64 


0.55 


0.59 


0.64 





0.3 


0.5 


0.7 


1 




0.63 


0.62 


0.64 


0, 


.66 


0.52 


0.57 


0.65 





0.17 


0.33 


0.67 


1 




0.54 


0.51 


0.55 


0, 


,65 


0.69 


0.52 


0.55 





0.05 


0.15 


0.5 


1 




0.35 


0.38 


0.34 


0, 


,53 


0.59 


0.30 


0.36 





0.125 


0.25 


0.5 


1 




0.33 


0.35 


0.31 


0, 


,67 


0.76 


0.27 


0.34 





0.1 


0.2 


0.5 


1 




0.34 


0.36 


0.33 


0, 


,61 


0.72 


0.29 


0.35 





0.3 


1 








0.05 


0.06 


0.05 


0, 


,75 


0.67 


0.04 


0.06 





0.4 


1 








0.19 


0.23 


0.18 


0, 


,72 


0.51 


0.15 


0.21 





0.5 


1 








0.45 


0.50 


0.42 


0, 


,58 


0.34 


0.36 


0.47 





0.6 


1 








0.71 


0.73 


0.70 


0, 


,39 


0.20 


0.63 


0.72 





0.7 


1 








0.78 


0.72 


0.79 


0, 


,22 


0.11 


0.78 


0.76 



aim is to estimate the malformation rate and tertiary aim is to estimate the overall toxic- 
ity rate as accurately as possible. To incorporate the multiple objectives in the study, one 
may foUovsf the strategy laid out in Cook and Wong (1994) to find an optimal design that 
provides user-specified efficiency for each objective. Clearly, the optimal design sought 
should provide higher efficiencies for more important objectives and user-specified effi- 
ciencies reasonable enough so that the optimal design exists. For space consideration, we 
do not provide multiple-objective optimal designs for simultaneously estimating the ef- 
fective dose for prenatal death rate, malformation rate and overall toxicity rate, but note 
that the key idea for finding such a design is to first formulate each objective as a convex 
function of the design information matrix and then combine all the convex objectives into 
a single convex functional using a convex combination. As described in Cook and Wong 
(1994), each set of weights used in the convex combination can be judiciously chosen to 
satisfy the efficiency requirement for each objective. In the case of a two-objective design 
problem, the weights and the dual-objective optimal design can be determined graphi- 
cally via efficiency plots. Wong (1999) provided several illustrative applications of such 
ideas to construct multiple-objective optimal designs in several biomedical problems. 

One may be rightly concerned that the optimal designs are dependent on the para- 
metric models. This dependence is inescapable but as we have advocated all along, the 
user must check robustness properties of optimal designs to all assumptions before the 
design is implemented. We focused on the simpler situation when we were concerned 
about misspecification of initial values, but if there is concern about other aspects in the 
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model assumptions, a similar strategy can be applied. For instance, one may question the 
validity of the WeibuU models to describe the malformation and prenatal death rates. 
If scientific opinion suggests alternative models may be more appropriate, one can then 
construct optimal designs for different models and compare their efficiencies under the 
competing models. The hope is that there is a design that remains efficient under all 
models on which experts may agree. 

Here is a short illustration of the situation just discussed: Assume, as before, that both 
the malformation and prenatal death rates have the same form and can be described using 
two plausible models: 

7^2 ('^) = 1 



1 + 62 rf"^^ 



and 



(3)^^^_1 «2 



-62+72^ ' 



Suppose the sets of initial values are 6*^ ' = (0.88,0.25,2.8), 6*^^^ = (0.91,4.3,3.5), 6*^ ^ = 
(0.94,1.3,5.1) and 6*^^ = (0.98,3.5,3.2). We recah that 6*^^^ = (0.06,0.7,3.37) and 9^^^ = 
(0.13, 0.27, 3.33). Here the superscripts denote the three different models used to describe 
the probability rates. 

Table 13 lists the locally optimal designs for a = 0.05 and their efficiencies under dif- 
ferent assumptions on the probability models. The robustness properties of each optimal 
design under each set of probability models can be compared. For this setup, the effi- 
ciency results are quite reassuring because the smallest efficiency in the table is at least 
0.76. Of course, different assumptions on the sets of initial values may not yield the same 
conclusions. 

In summary, our proposed design strategy is quite general and possess several advan- 
tages over existing methods. Unlike uniform designs, our approach is based firmly on 
statistical principles and the proposed maximin optimal design provides good protection 
against misspecification in the initial values of the model parameters. The optimal design 
allows prior information to be included in its construction and, if required, can also incor- 
porate multiple objectives with possibly unequal interests. Consequently, the proposed 
optimal design is able to meet the practical needs of the researcher more adequately than 
current designs. 

Appendix: Proofs 

Proof of Theorem 1. From (3.1), the information matrix for a design ^ can be repre- 
sented as 

k 

M{^,e)^Y.'^j{d.,)f^{d.o, 

i=l 
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Table 13. Various locally optimal designs (leftpart) and their efficiencies under different prob- 
ability models for estimating of prenatal death (first three rows), malformation rate (middle 
three rows) and overall toxicity (last three rows) 



Model 




X2 


X3 




W2 


Wo, 


^(1) 


^(2) 


^(3) 


weibull 





0.686 


1 


0.396 


0.548 


0.056 


1.000 


0.910 


0.869 


model 2 





0.624 


1 


0.417 


0.546 


0.037 


0.911 


1.000 


0.968 


model 3 





0.630 


1 


0.354 


0.561 


0.086 


0.853 


0.940 


1.000 


weibull 





0.616 


1 


0.297 


0.602 


0.101 


1.000 


0.954 


0.832 


model 2 





0.662 


1 


0.284 


0.592 


0.123 


0.918 


1.000 


0.503 


model 3 





0.535 


1 


0.290 


0.610 


0.100 


0.882 


0.768 


1.000 


weibull 





0.654 


1 


0.323 


0.550 


0.127 


1.000 


0.885 


0.726 


model 2 





0.581 


1 


0.357 


0.521 


0.121 


0.825 


1.000 


0.910 


model 3 





0.544 


1 


0.297 


0.564 


0.138 


0.761 


0.942 


1.000 



where the vector / is defined by 

D{d) 



fid) 



Vm(l + (m - l)(j){d))n2{d)il - 7r2(d)) 



i^-^2{d)) 



1 



m(rf)(l + (m(d) - l)0(d))7r2(d) W^d''^ \n{d) j ' 



We now apply Elfving's theorem (see Elfving (1952)), which gives a geometric charac- 
terization of the optimal design. More precisely, from this result it follows that a design 
^ = {di,Wi}^^i is locally optimal if and only if there exist numbers ei,. . .,ek S {—1, 1} 
such that for some i/ G K the point 

z.P = ^(^0,l/62,^ln(^-i^^ii-^))'' = ^e,z«,M) (A.l) 

is a boundary point of the Elfving set 

7^ = conv({e/(d) | d G [0, 1], e G {-1, 1}}). (A.2) 

A typical picture of this set is presented in Figure A.l for the case of constant functions 
(f) and m. We note that the curve 

A' = {/(d),dG[0,l]} 
is contained in subspace {x = {xi,X2,X3)'^ G R'^| xi > 0} and the set 

{(l,dT%62d^^ ln(d))|dG [0,1]} 
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-0.04 



-0.064 



-0.02 



0.02 



0.04 



0.061 







-2 




-0.5 



-1 



Figure A.l. The Elfving set defined in (A. 2) for tlie parameters a = 0.133, b = 0.272, 7 = 3.33. 
Tiie points /(di), —/(di) and /(da) are denoted by A, C and B, respectively, while the point 
i^P is defined in (A.l). 



defines a U-shaped curve. From the monotonicity assumption for the function (3.3), it 
follows that the curve X is also U-shaped (see also Figure A.l). We denote the end-points 
of this curve by A and B and recall that the first coordinate of the vector P is equal to 
and that f is the scaling constant such that I'P touches the boundary of the Elfving 
set TZ. Note that in the case a — > we have that 



for some constant c and, consequently, the vector vP touches the boundary at the plane 
£ spanned by the points A, B and C, where A and B correspond to the doses and 
1, respectively, and — C corresponds to a third dose, say d* € (0,1). Consequently, the 
locally optimal design is a three-point design with support points 0, 1 and c?*, if a is 
sufficiently small. In the case where a — s- 1, the situation is exactly the same and the 
locally optimal design is also supported at three points, including the boundary points. 
From the geometry of the Elfving set TZ we see that there is also direction P, where the 
intersection with the Elfving set can be represented by two points of the curves X and 
—X. In particular, this situation occurs if a « 1 — c"''. In this case we have 
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for some constant c and the locally optimal design is supported by a two-point design. 
Moreover, if a moves from to 1, it follows from the geometry of the Elfving set that 
the situation is changing continuously, which proves the assertion of the theorem. □ 

Proof of Theorem 2. Multiplying the last coordinate of the vector equation (A.l) by 
72 yields 



vP 




to(1 + {m - l)(j))Ti2{di) 




(A.3) 



for the boundary point vP G TZ. If {(i*(l); w*} denotes an optimal design for the param- 
eter 9 = (02,62, 1), it follows that equation (A.3) holds for this design with 72 = 1. Now 
it is easy to see that (A.3) is also true for the design {(d*(l))-^/'''2; w*} for the parameter 
6 = (02,62,72), where 72 > is arbitrary. This proves the first part of the theorem. The 
second part follows by similar arguments, which arc omitted by the sake of brevity. □ 

Proof of Theorem 3. Part (a) of the Theorem follows directly from the geometry of the 
Elfving set. If the locally optimal design is supported at three points, the corresponding 
point vP touches the Elfving set in the plane spanned by the points A, B and C, which 
does not depend on the value of a. For a proof of part (b) we note that, according to 
Elfvings theorem, a locally optimal design of the form {0,(^2, wi,W2} must satisfy the 
equation 



In 





1/62 

ln(l - a) 



ewig{Q) I 



ZW2g{d2) 



1 

b2df ln(d]=) 



where the function g is defined by 



9{d) 



(l-^2(d)) 



m{d){l + {m{d)-l)(j){d))^2{d) 
It is easy to see that this equation yields 



1/62 



V I In(l-a) 



= -£W2\ 



(l-7r2(d2)) 



m(d2)(l + (m(d2) - I)0(d2))7r2(d2) V^sd? Hdf) 



df 
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which simphfies to the equation 



1 \ _ / (1-712(^2)) / 1 

HEDZ^)) '""'V mid2)il + im{d2)-md2))7:2{d2r ' {HdT) 



It foUows that d2 ~ ED a. Since wi ^ \ — W2 from the equahty for the first coordinate, 
we have that W2 = g{0)/{g{0) + g{EDa)). □ 
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